On the Q-Ball Profile Function 
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We use analytic and numerical methods to obtain the solution of the Q-ball equation of motion. 
In particular, we show that the profile function of the three-dimensional Q-ball can be accurately 
approximated by the symmetrized Woods-Saxon distribution. 

1 Introduction 



The existence of Q-baUs is a general feature of scalar field theories carrying a conserved 
U{1) charge [1, 2]. Q-balls can be understood as bound states of scalar particles and appear 
^ . as stable classical solutions (nontopological solitons) carrying a rotating time dependent 
internal phase. They are characterized by a conserved nontopological charge Q (Noether 
charge) which is responsible for their stability (see, for example, Refs. [3]- [4]). These 
CN ■ features differentiate the Q-ball interaction properties from those of the topological solitons 
since here the charge Q can take arbitrary values in a specific range, allowing for the 
I possibility of charge transfer between solitons during the interaction process. 

Up till now, comprehensive studies of these objects have been made by using either 
numerical simulations [3]- [5] or some analytic considerations [1, 6, 7]. Recently, in [8] the 
Qh! explicit relation between the energy and the charge of the Q-balls has been derived using 
r-| ! analytic arguments. In particular, a semi-Bogomolny argument in the energy density led 
to a first order differential equation whose solution however, did not satisfy the correct 
boundary conditions and differed considerably from its exact form. In this work, we 
present a method to obtain an analytic form for the Q-ball profile function which is in 
I good agreement with the numerical results. 

We consider the U{1) Goldstone model, given by the Lagrangian 

c = ^d^(bd^^-um) (1) 

where is a single complex scalar field in three spatial dimensions while the potential U{\(j)\) 
is a function of \(f>\ only and has a single minimum at = 0. This is equivalent of stating 
that there is a sector of scalar particles (mesons) which carry t7(l) charge and have mass 
squared equal to ^U"{0). The corresponding energy functional is given by 

E^J Q|a,<^p + i|V0p + UM)^ dPx. (2) 

The model has a global U{l) symmetry and an associated conserved Noether current 
exists whose covariant conservation 9^ = leads to the existence of the conserved Noether 
charge Q given by 

Q^Yil {^9tcl>^cl)dt^)d^x. (3) 



^Permanent Address: Institute of Mathematics, University of Kent, Canterbury CT2 7NF, UK 



A stationary Q-ball solution has the form 



e 



fir) (4) 



where /(r) is a real radial profile function which satisfies the ordinary differential equation 

f"{r) + lnr) = -u;y{r) + U'{f) (5) 

with the condition /(oo) = and /'(O) = 0. This equation can either be interpreted as 
describing the motion of a point particle moving in a potential with friction [1], or in 
terms of Euclidean bounce solutions [9]. In each case the effective potential being Ueff(f) = 
ijp-p 12 - leads to constraints on the potential f/(/) and the frequency w in order for a 
Q-ball solution to exist. Firstly, the effective mass of / must be negative. If we consider a 
potential [/(/) which is non-negative and satisfies t/(0) = t/'(0) = 0, ?7"(0) = > then one 
can deduce that w < Furthermore, the minimum of U{f)/p must be attained at some 
positive value of /, say < /o < oo and existence of the solution requires that uj> uj- where 
= 2C/(/o)//o. Hence, Q-balls exist for all w in the range < < 

Then, the charge and the energy of a stationary Q-ball solution simplify to 



Q = Attu 



Jr''f{r)dr (6) 

E = 47ry Qa;V'W + ^/"(0+f/(/))r^rfr. (7) 

It has been observed using numerical and analytic methods that the classical stability of a 
Q-ball is related with the dependence of its charge on the internal frequency lj. For small 
internal frequency, close to its minimal value w_, the profile function is almost constant 
which implies that the charge (6) is large and this corresponds to the so-called thin-wall 
approximation. On the other hand, for large internal frequency (close to its maximal value 
u+) the profile function (and thus the charge) tends to zero and this corresponds to the 
thick- wall approximation. In fact, for lo ^ u!+, the behavior of the charge Q depends on 
the particular form of the potential and the number of dimensions [7] . In the case studied 
here we show that Q ^ oo as u ^ uj+. 

The choice of the potential is not unique; the standard requirement is that the 
function U{f)/p has a local minimum at some value of / different from zero. For simplicity, 
we consider the following form 

U{f) = /2 (1 + (1 - ff) (8) 

which implies that in terms of the earlier notation we have that w+ = 2 and w- = V2, and 
therefore stable Q-balls exist for \/2 < uj < 2. In this case, the energy-charge dependence 
of the Q-balls, obtained in [8] using a semi-Bogomolny argument, is given by the analytic 
expression 

Enog - y'^QBog + 27/6 ^ Bog + ^1176^ ^ Bog + 2592 2V6 3l/3 ^ Bog 



1,^1^ (20 - 547r^ + 277r4) 
1944 25/6 32/3 



QLf + (9) 



For large Q, the above equation corresponds to the upper energy bound. The same result 
can be obtained by representing the Q-ball profile function by the Woods-Saxon distribu- 
tion (ie a generalization of the semi-Bogomolny solution) as presented in [8]. 



Let us emphasize that neither the semi-Bogomolny nor the Woods-Saxon ansatz 
describe accurately the Q-ball profile function. In particular, in both cases the derivative of 
the profile function at the origin is non-zero (that is, different from the required boundary 
conditions) while the corresponding profile function differs from the exact one obtained 
by solving numerically the equation of motion. In the next section, we show that the 
symmetrized Woods-Saxon distribution describes accurately the Q-ball profile function 
and satisfies the correct boundary conditions. 



2 Equation of Motion 

The Q-ball equation of motion (5) for the specific potential (8) is given by 

fir) + \ fir) = (4 - u') f{r) - 8 fir) + 6 fir) (10) 

and satisfies the boundary conditions /'(O) = 0, /(oo) = 0. For large values of the argument 
r the nonlinear terms in (10) can be neglected and we get the asymptotic behaviour /(r) ~ 
exp(-\/4 - iv'^r). The exact solution of (10) should contain only one free parameter iv so 
that the charge Q the energy E and the initial value /(O) should in principle be expressed 
as functions of this parameter. Since equation (10) is too complicated to be tackled by 
analytic methods, we shall seek a suitable ansatz for the profile function and search for 
relations connecting E,Q and /(O) with uj. 

As mentioned earlier, a test profile of the Woods-Saxon type cannot reasonably 
approximate the Q-ball solution since it fails to satisfy the correct boundary conditions. 
Nevertheless, it implies an energy-charge relation (9) which holds remarkably well for a 
wide range of energies, a fact that certainly requires some further investigation. A more 
general test function that satisfies the correct boundary conditions and has the correct 
asymptotic behaviour is given by the symmetrized Woods-Saxon distribution 

= \f , • 

+ Ci cosh (ar) 

The values of the arbitrary parameters c, ci and a can then be determined by fitting the 
data of the numerically solved (10). Having done that, we got a very satisfactory agreement 
in all cases. It is interesting to realize that /(r) satisfies the following differential equation 

/"(-) + ;/v) = T (' - i) - S (' - i) + 1^ + 3I) + ° <^') <'^' 

which in the limit ar > 4 looks exactly like (10). It is then quite reasonable to expect 
that there must be a critical value of r = rc beyond which the following equations must 
approximately hold true: 

(13) 



4 \ ar, 

The thin-wall approximation is reached when ci ^ so we will assume that ci < 1 and we 
shall neglect the last term in the left hand side of (15). Then the system above gives the 



approximate solutions 



(16) 



4V2 



3 4-^ 



(17) 
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_ 1 , /1 3 4- 
2 + Y ^ 



(18) 



A comparison with values obtained by a direct fit shows that the relations above are quite 
good. In fact, equation (16) is excellent, while equation (17) gives the right shape but the 
actual values for a are on the average 10% higher than expected, implying a faster drop of 
the profile function. Note that both c and a are slowly varying functions of w. 
The equation of motion (10) can be written as: 



(19) 



In the absence of the friction term, equation (19) would simply imply the conservation of 
energy for the corresponding mechanical problem. In the presence of friction and upon 
integrating (19) we get that the initial potential energy equals the work done by friction. 
This relation can be used to provide a further constraint on the form of (11) 



-a;V'(0)-C/(/(0)) 



Jo 



r{r) 



dr ■ 



(20) 



Recall that /(O) = c/^/^ + Cl. Finally, using the symmetrized Woods-Saxon distribution (11), 
the charge (6) and the energy (7) of the Q-ball can be explicitly evaluated in terms of the 
parameters c, ci and a: 



Q = imo 



cosh 



TT + cosh 



E 



An 



a2+4(4 + a;2) / \ / a' , „A ' 



where the integrals /q, h and h are given by 

Q 



h = 



2 d, 
C Ci—Io, 

aci 



c% + 2cV— Jo 
aci 



1 



(21) 
(22) 

(23) 
(24) 

(25) 



It is readily seen that the actual values for Q and E depend strongly on the value of ci 
especially for small values of ci. In this region, ci can be eliminated in favour of Q to get 

(I 24 + 3^2 1 , 24 + 3a2 ,^1 id,+ a'^)'n:'^ 

E = Q + ^6^M-Q^ + ^^M-Q^- ^ 8a; + ^iQ ^) (26) 



where 



It is interesting to see that the leading terms of (26) are identical to those of (9) for a = 2V2, 
c = 1, u = \f2. Taking into account the fact that a and c are slowly varying functions of w 
helps to explain the wide range of applicability of (9). 

As a result, the parameters c, ci, a and uj for specific values of the charge Q can 
be uniquely determined by solving the system of equations: (16), (20), (21), and (22) 
provided that the energy-charge dependence of the Q-ball is given by (9). We expect the 
values obtained that way to be quite accurate as long as (9) remains reasonably accurate. 
In practice, this means that Q > 16, a value far distant from the values normally associated 
with the region of the thin-wall approximation. In Table 1 the values of the arbitrary 
constants in (11) are presented for different values of Q (or w). That way, the initial value 
/(O) can be determined explicitly; a comparison with values obtained numerically show 
that the symmetrized Woods-Saxon distribution describes accurately the Q-ball profile 
function. 
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c 
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/(O) 


f {^\iu7n 


23.68 


1.79 


1.097 


0.12 


2.66 


1.037 


1.024 


35 


1.72 


1.087 


0.068 


2.60 


1.052 


1.055 


61.6 


1.64 


1.068 


0.0243 


2.54 


1.056 


1.065 


149.3 


1.57 


1.049 


0.0031 


2.53 


1.048 


1.056 



Table 1: Values of the parameters c, a, a and /(O) obtained from (11) and numeri- 
cally (/(0)„„to) for different values of ui. 

Finally, Figure 1 presents the profile function obtained analytically and numerically 
for ui = 1.64. Figure 2 presents the values of the profile function at the origin for different 
w obtained numerically and analytically. 

3 Conclusions 

In this work, the basic properties of the Q-ball profile function have been extensively 
studied by means of mainly analytic methods. In particular, it has been shown that 
the profile function can be accurately approximated by the symmetrized Woods-Saxon 
distribution, while the corresponding energy and charge can be explicitly evaluated. The 
approach presented here might prove to be particularly useful in understanding the basic 
properties of Q-balls such as existence, small vibrations and stability. We believe that a 
similar line of argument can be applied to study the profile function and the energy-charge 
dependence in other types of potentials. In fact, we expect that the symmetrized Woods- 
Saxon distribution will accurately describe all types of Q-balls independently of the exact 
form of the scalar potential. Work in this direction is currently in progress. 
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Figure 1: The profile function /(r) as a function of r for uj = 1.64 




Figure 2: The frequency dependence of the profile function at the origin /(O). 
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